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We review the experimental searches for new particles in the dijet mass spectrum con- 
ducted at the CERN SppS, the Fermilah Tevatron Collider, and the CERN Large Hadron 
Collider. The theory of the QCD background and new particle signals is reviewed, with 
emphasis on the choices made by the experiments to model the background and signal. 
The experimental techniques, data, and results of dijet resonance searches at hadron 
colliders over the last quarter century are described and compared. Model independent 
and model specific limits on new particles decaying to dijets are reviewed, and a de- 
tailed comparison is made of the recently published limits from the ATLAS and CMS 
experiments. 

1. Introduction 

Experiments at hadron colliders have used the dijet mass spectrum to search for new 
particles beyond the standard model. At the CERN SppS, the Fermilah Tevatron 
Collider, and the CERN Large Hadron Collider, with each successive advance in 
collision energy and integrated luminosity, progressively more energetic collisions of 
the partons in the incoming hadrons are produced and observed. Each machine in 
its time has therefore probed the highest masses of dijet resonances: new particles 
that decay into two partons, giving two jets in the final state. The simple process 
the experiments have searched for is the s — channel production and decay of dijet 
resonances shown in Fig. [1] 

Here we review these experimental searches, their techniques, data, results, and 
limits on new particles. In section [2] we review the theory of the QCD background 
and the models of new particle signals. In section[3]we review the experiments from 
each collider in chronological order. In section [3.11 we review how each experiment 



*Also affiliated with Fermilab during this review. 



2 Robert M. Harris and Konstantinos Kousouris 




Fig. 1. Diagram of dijet resonance in the s — channel. The initial state and final states contain 
two partons (quarks, anti-quarks or gluons) and the intermediate state contains a resonance X. 

modeled the resonance shapes as a function of the dijet mass. In section 13.21 we 
review the data of each experiment and how each experiment modeled the QCD 
background. In section [3731 we review the limits on dijet resonance masses published 
by each experiment, discussing the experimental uncertainties, statistical proce- 
dures, and the cross section assumed for each model. In section 2] we conclude with 
a few observations. Also, in [Appendix A[ we include details of the cross-section cal- 
culations for axigluons and excited quarks by ATLAS and CMS, which are necessary 
to understand the mass limits on these models from the two experiments. 

2. Theory 

In this section we present the fundamental ingredients of the theory, which are 
necessary for the better understanding of the experimental searches presented in this 
review. In Section[5TT]we describe some basic features of Quantum Chromodynamics 
(QCD), and in Section [52] we present the theoretical models that predict partonic 
resonances and are quoted in the experimental searches. It should be noted, that 
the purpose of this section is not to give all the details of the models presented, but 
rather an overview of their features. 

2.1. Elements of QCD 
2.1.1. The QCD Lagrangian 

Quantum Chromodynamics is the gauge field theory of the strong interaction be- 
tween particles that carry the color degree of freedom. The underlying symmetry 
group is the SU(3)c which makes QCD a non-Abelian theory. The profound im- 
plication of this property of QCD is that the gauge mediators (gluons) are colored 
and thus self interacting. The QCD Lagrangian is written as: 

CqCD = q^,a (ird^Sab ^ .9.7^*^,0^1 " m.Sab) q^^ ~ -^F^uF''"'^ , (1) 

i 

where ^ represents the quark spinor of flavor i and color a = 1 — > 3, G^^ is the 
gluon field associated with the generator (A = 1 H> 8), is the gauge coupling, 
and F^^ is the gluon field tensor: 

F^v = d^iG^ - dyC^ - gsfABcG^C'^ . (2) 
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The structure constants Jabc satisfy the relation: 

[t^,t^]=^.fABct'^■ (3) 

The non-Abehan nature of QCD leads to two remarkable features: the confine- 
ment and the asymptotic freedom. As a result of the confinement, only color-singlct 
states can be directly observed, which means that quarks and gluons cannot be 
found free. The asymptotic freedom is the property where the running strong cou- 
pling constant decreases with increasing momentum transfer between the strongly 
interacting particles. This in turn means, that the hard-scattering of quarks and 
gluons can be described in a perturbative way. 

2.1.2. Formation of jets 

Because of the confinement, partons cannot be detected free. Instead, the experi- 
mental signatures of quarks and gluons are the jets. A jet is a "spray" of highly 
coUimated particles, primarily hadrons, but also photons and leptons. A jot is a not 
uniquely defined object, but the output of a well-defined mathematical rule (clus- 
tering algorithm) , which clusters the jet constituents, according to their kinematic 
properties. This procedure is based upon the features of QCD, which describe the 
transformation of a parton to a set of observable particles. The jet-formation steps 
are the following: 

• Parton branching: each parton, whether a gluon or a quark, has a fi- 
nite probability to split into two partons, which are emitted in small an- 
gles with respect to the direction of the initial parton. One feature of 
the parton branching is that the probability depends on the color factor 
related to the type of the involved partons. For gluon — > gluon, gluon, 
gluon — > quark, antiquark, and quark — >■ gluon, quark splittings, the color 
factors are Cgg = 3, Cgq = |, and Cqg = ^ respectively. As a result, glu- 
ons systematically shower more than quarks. Another implication of parton 
branching in small angles is that throughout the process, partons are pro- 
duced close to the direction of the initial partons, which results in a high 
degree of collimation of the final hadrons. It should be noted, that the par- 
ton branching is a perturbative procedure, which can be re-summed to all 
orders of the perturbation series, under certain assumptions. 

• Hadronization: when the parton shower has evolved long enough, the 
energy of the partons is reduced, such that low-momentum transfer oc- 
cur. In these conditions, the parton interactions become non-perturbative, 
and the phase of hadronization begins. During the hadronization, par- 
tons are combined into color-singlet states, thus forming the hadrons. 
While the hadronization procedure cannot be described perturbatively, the 
local parton-hadron duality ensures that the flow of quantum numbers at 
the hadron level, follows approximately the corresponding flow at parton 
level. 
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• Underlying event & out- of- cone showering: the term underlying event 
in hadron collisions is used to describe the activity not related to the hard 
scattering, for example due to multiple parton interactions happening si- 
multaneously. Since the definition of jets involve the clustering of hadrons 
which are sufficiently correlated, it can happen that particles originating 
from the soft interactions are clustered together with those coming from 
the hard-scattered parton shower. In the opposite direction, partons from 
the initial shower can be emitted in relatively large angles, and the asso- 
ciated hadrons may not be clustered in the resulting jets. This effect is 
commonly known as out- of- cone showering. 

Despite the fact that the formation of jets is a complicated effect, certain con- 
clusions can be drawn: to first approximation, the kinematical properties of a jet 
are the same as those of the original parton. However, the various effects involved, 
introduce an intrinsic resolution of the hadronic jet properties with respect to the 
parton properties. 



2.1.3. Kinematics of two-parton scattering 

Before the details of the strong interaction dynamics are discussed, it is useful to 
present the kinematical properties of a two-to-two parton scattering. In the topology 
of the 1 + 2^-3 + 4 scattering, some general kinematic relations hold, regardless of 
the details of the interaction. The Mandelstam variables of the process are defined 
as s = (pi -\-p2)'^i t = {pi — Ps)^, and u = {p2 — Ps)^, where pi are the four-momenta 
of the partons. For massless partons, the Mandelstam variables satisfy the relation 
s + i + M = and two of those can be expressed as a function of the third one and 
the scattering angle 9* in the center-of-mass frame: 

f= -is(l-cosr), u = -is(l + cosr) . (4) 

The rapiditie^ of the outgoing partons, in the center-of-mass frame, are opposite 
(±y*), due to transverse momentum conservation, and related to the scattering 
angle: 

cos 6** — tanhy*. (5) 

The Mandelstam variable s can be expressed in terms of the outgoing partons 
transverse momentum px and y*: 

s = 4p|cosh2y*. (6) 

In the laboratory frame, the rapidities 1/3.4 of the outgoing partons are related to 
the rapidity of the center-of-mass frame y and to y* as: 

y= — — > y = — — • (7) 
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From the relations above, one can express the scattering angle at the center-of-mass 
frame as a function of the rapidities of the scattered partons at the laboratory frame: 

cos6l* tanh (8) 



2 

Each initial parton is carrying a fraction x of the hadron momentum and the in- 
variant mass of the two-parton system is expressed as: 

— s — X1X2S, (9) 

where Xi^2 are the momentum fractions of the interacting partons and is the 
colliding energy of the hadrons: 

xi = coshy*e^, X2 = coshy*e^^. (10) 

Following from the relation above, the rapidity of the center-of-mass frame y can 
be expressed as a function of the momentum fractions: 

2/= Jin- (11) 

2 X2 

2.1.4. Partonic cross sections 

The dynamics of the hard scatter of colliding hadrons is approximately described as 
a two-to-two process between massless partons. Because of the different structure 
and color factors of the interaction between the parton types, the matrix elements 
are different for each subprocess. The leading order (LO) amplitudes can be calcu- 
lated analytically using the Feynman diagrams at tree level, and are summarized 
in Tabic [TJ The squared amplitudes are averaged (summed) over the initial (final) 
color and spin indices, and are expressed in terms of the Mandelstam variables. 

Figure [2] shows the matrix elements of the various subprocesses, at LO, as a 
function of cos 6**. With the exception of one subprocess {qiqi -> 9292), there is a 
characteristic t — channel pole which enhances the two-parton scattering at small 
angles. Another important feature is the fact that, due to the larger color factor of 
gluons, the matrix element of the subprocesses with gluons in the initial state lead 
to larger values. 

2.1.5. Hadronic cross sections 

In a hard scattering process, initiated by colliding hadrons, the experimentally mea- 
sured cross section can be generally expressed in terms of the parton distribution 
functions (PDFs) f{x) and the parton-parton scattering cross section ct, summed 
over all the incoming and outgoing parton types (because the experimentally ob- 
served jets cannot distinguish between the parton types) : 

0- = ^ / dxidx2fiixi, iJ,%)fj(x2, Hp)aij asinj^), — , — ) (12) 
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Table 1. Leading order matrix elements for 
2-to-2 massless parton subprocesses from Ref. [T] 
The color and spin states are averaged over the 
initial states and summed over the final ones. 



Subprocess 



9192 -S- 9192 
9192 9192 

99 99 
91 91 -S> 9292 
99 99 
99 -s- 99 
99 99 
91 91 
99 -s- 99 



4 J^+u^ 



9 P 
4 s^+u^ 



9 £2 



9 V t2 ^ ) 27 lit 



9 ~~g^ 



/ 27 It 



32 P+v? _ 8 P+v? 
27 tfl 3 3^ 

1 t"2+fl2 _ 3 t^+uf. 
6 t" u 8 ^2 

4 s^+ii^ , fl^+s^ 



9 su 



9 / q tli su st 

2 V fl' 



In the equation above, Q is the characteristic hard scale of the interaction (e.g. the 
dijet invariant mass in a two-to-two parton scattering), is the factorization scale, 
which is of the same order as Q and separates the long-distance, non-perturbative 
interactions from the hard scattering, and ^ib. is the renormalization scale. Both 
the ^F,R scales are arbitrary parameters of a fixed-order calculation. At all orders 
of the perturbative expansion, the cross section should be independent of them 
[da/diiR = da/d^F = 0). In all practical calculations of cross sections at a fixed 
order, it is assumed that fiji — fip ^ Q. It should be noted, that the higher the 
order of the calculation, the weaker is the dependence on (ir.f- 

It is often helpful in hadron collisions to quantify the effect of the parton distri- 
bution functions by introducing the parton luminosity factor. This is defined as: 

= J J dxidx2fi{xi)fj{x2)S{xiX2 - t), (13) 

where 

T = X1X2 = -. (14) 
s 

In practice, experimental constraints are imposed on the rapidities of the out- 
going partons, observed as hadronic jets. It is therefore more convenient to express 
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cose* 



Fig. 2. Leading order matrix elements for two-to-two massless parton scattering, as a function of 
cos 9* . 



the parton luminosity as a functions of the variables r and y, rather than 2:1^2: 
■drdy = 

dLij {ymin ; yvaax) 



dxidx2 = ^^rf^drdy = drdy. The parton luminosity then is: 



dr 



h (V^e^) {V^e-y) dy 



(15) 



The hadronic cross section of any process, can be expressed generally as a function 
of the parton luminosity factor and the partonic cross section!^ 





1 dLij 




s dr 



(16) 



In the specific case of a two-to-two scattering, resulting in the production of 
two jets, the differential cross section can be expressed as a function of the di- 
parton invariant mass and the scattering angle at the center-of-mass frame. To 
first approximation, this cross section is equal to the observed dijet cross section. 
The matrix elements presented in table [1] are folded with the parton distribution 
functions, giving: 



d'^ahad Tra^ 



dm d cos 6* 



E 



1 dL,j 
s dr 



T—ni^ / s 



Tij {cos 0*), 



with 



J"ij(cosr) = "^S{ij kl) — ^- 

kl ^ 



Ski 



(17) 



(18) 



In the equations above m \/§ — ^/rs is the dijet invariant mass. 
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2.2. Models of partonic resonances 

2.2.1. Excited quarks 

In various theoretical models, ordinary quarks can be composite objects with a 
characteristic compositeness scale A. As a natural consequence, excited states are 
expected, called simply excited quarks and denoted by q* . Depending on the details 
of the composite models, the excited quarks can have various values of spin and 
weak isospin. In the simplest case, they take the value of 1/2. The interaction of 
excited quarks with the Standard-Model fields is of a "magnetic" type, and the 
Lagrangian takes the form SI 

C = ^rR'j'-'' (^gsfstaG'^^, + gf-^ + 9'r\B,}j qL + h.c, (19) 

where ta and r are the generators of the color SU{3) and isospin SU{2), Y is the 
hypercharge, G"^, W^^^, B^i, are the field tensors, gs,g,g' are the gauge couplings, 
and /s,/, /' are dimensionless constants, accounting for possible deviations from 
the Standard-Model couplings. 

In hadron collisions, the production of an excited quark happens through the 
quark-gluon fusion. Subsequently, the excited quark decays to an ordinary quark 
and a gauge boson. The dominant decay channel is q* — ^ qg, leading to a dijet 
signature. The partial width for the decay of an excited quark with mass m* is 
given by: 

r{q*^qg)^^asf^j^ (20) 

2.2.2. Randall- Sundrum gravitons 

The gravity model from Randal and Sundrum (RS) was proposed as a solution 
to the electroweak vs Planck scale hierarchy problem. In this model the hierarchy 
is generated by an exponential function of the compactification radius of one extra 
dimension. The metric in the 5-dimensional space is given by: 



= e-2fe'^c0^^^2;^2.i^ ^ ^2^02^ ^21) 

where (p is the extra dimension with compactification radius Tc, fc is a constant of 
the same order and related to the 5-dimensional Planck scale M, and are the 
usual space-time dimensions. The reduced effective 4-D Planck scale Mpi is given 
by: 

Ml, = ^ (1 - e-^'^-^) . (22) 

In this model, spin-2 gravitons appear as the Kaluza-Klein (KK) excitations of 
the gravitational field h^^ , whose coupling to the Standard-Model fields is given by 
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the interaction Lagrangian H 

1 

A 



>C/ - - — h'^'-T^,, (23) 



with T^i^ being the energy-momentum tensor of the matter fields. The scale and 
the mass m„ of the KK excitations can be expressed as a function of the fundamental 
constants of the theory: 

= Mpie-""^^^ , m„ = fcx„e-'='^=". (24) 

The coupling constant of the graviton-matter interaction is the inverse of the scale 
A^: 

A^ m„ 

where x„ is the n-th root of the the Bessel function of order 1. The phenomenological 
consequences of the RS-gravitons are essentially determined by their mass, and the 
ratio k/Mpi. If the fundamental constants of the model satisfy the relation krc ^12, 
then Att ~ TeV, and RS gravitons can be produced in hadron collisions. Through the 
graviton coupling to the matter fields, RS-gravitons can decay to partons, leading 

Q 

to a dijet signature. The relevant partial widths ° for the first KK excitation are 
given by: 



9 

xi 



2 



and 



T{G ^ qq) = ^ ( mi. (27) 
^ leOvr \MpiJ ^ ' 



2.2.3. Axigluons 

In the axigluon models the symmetry group of the strong sector is expanded 
to a chiral color group SU(3)l x SU(3)i{ which, at some energy, breaks to the di- 
agonal SU(3). Under such a symmetry group, the left-handed and right-handed 
fermions "01-, -R = ^(1 T 75)'0 transform differently and the transformations are gen- 
erated by the r£ ^ generators. Equivalently, the group can be described by a linear 
transformation of the generators, divided into vectorial Ty — + and axial 

= T£ — Tj!; ones. The associated gauge field to the vectorial generators is iden- 
tified as the usual color field of QCD, while the gauge field associates to the axial 
generators is called the axigluon field. While the exact implementation of the chiral 
color group is model dependent, there are two universal features: the existence of a 
massive color octet axigluon field (corresponding to the broken symmetry) , and the 
existence of new particles which are needed to cancel the triangular anomalies. 

Axigluons can decay to quark-antiquark pairs, which leads to a dijet experi- 
mental signature. Note that due to parity conservation, the axigluon cannot decay 
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to a gluon-gluon pair (all gluon-axigluon vertices must have an even number of 
axigluons). The axigluon decay to fermions is described by the Lagrangian: 



La = -igs I 'l^^^lf.taqj I A""", (28) 

where gs = -s/IttoI, A is the axigluon field, and are the usua l color group gener- 
ators. The width of the axigluon decay can be shown to be^I^ 

r. ^ (29) 

where N refers to the open decay channels, and Ma is the axigluon mass. The latter 
is a free parameter of the theory, determined by the chiral color symmetry breaking 
scale and the details of the model. 



2.2.4. Colorons 

Similar to the axial color models, their exist other possibilities to enrich the grou" 



structure of the strong sector. Such a model is the flavor-universal coloron 
where the gauge group is extended to SU(3)i x SU(3)2. The corresponding gauge 
couplings are denoted as ^1^2 • Additionally, the model includes a scalar boson $, 
which develops a non-zero vacuum expectation value and breaks spontaneously 
the symmetry of the two groups. The diagonal subgroup remains unbroken and is 
identified as the familiar color group of QCD. In the rotated phase of the physical 
gauge fields, the initial gauge bosons are mixed, forming an octet of massless gluons 
and an octet of massive colorons. The mass of the colorons is expressed as a function 
of the fundamental parameters 



Mr = 



l) /, (30) 



where 9 is the gauge boson mixing angle with cot 9 = and / is the vacuum 
expectation value of the scalar field. The Lagrangian of the interaction between the 
colorons field and the quarks is similar to QCD: 

£ ^ -gs coi9 q^l^taq}^ (31) 

The above interaction predicts the decay of colorons to quarks with kinematically 
allowed masses. It can be shown that the decay width is: 

Tc « coe 9Mc: (32) 
6 



where N is the number of quark flavors with mass less than Mc/2. 
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2.2.5. Color octet scalars 

In various theoretical models, bosonic states can arise from gluon-gluon fiasion. 
The color octet scalar model {Ss) is one example of exotic color resonances ^1^. The 
coupling of the color octet scalar field with gluons is expressed with the Lagrangian: 

C = g,jd'^'"^S^,G%G'='^^-', (33) 

where gs is the strong coupling constant, k is the scalar coupling, A is the char- 
acteristic scale of the interaction, d"'"^ are structure constants of the SU(3) group 
defined by the relation {i",^''} = ^S"-^ + tP^H", and ^s, G^,^ are the color octet 
scalar field, and the gluon field tensor, respectively. The width of the color octet 
scalar resonance is given by: 



2.2.6. Z' and W 

New gauge bosons arise in models where the symmetry SU(2)i ® U(l)y of the 
electroweak Standard-Model sector is enlarged. Common features in these models 
are the new gauge coupling constants, which are of the same order as the SU(2)i 
coupling of the Standard Model, and the existence of new gauge bosons, namely 
W and Z' . Under the assumption that the new gauge bosons couple to ordinary 
quarks and leptons similar to their Standard-Model counterparts, the cross sections 
of these particles are calculated by scaling the the corresponding Standard-Model 
cross section. In particular, the Fermi constant Gp becomes 

where M and M' are the masses of W or Z , and W or Z', respectively. 



2.2.7. String resonances 

According to the string theory, particles are created by vibrations of relativistic 
strings, with mass Af^, and they populate Regge trajectories, which relate their 
spins and masses. In principle, the mass of the fundamental strings is of the order 
of the Planck scale. However, in some theories with large extra dimensions, it is 
plausible that Mg lies in the TeV scale. In this case, Regge excited states of quarks 
and gluons occur in hadron collisions. If the string coupling is small, the basic 
properties of the Regge excitations (production cross section and width) are model 
independent (the details of the compactification are irrelevant). This statement is 
true for parton scattering involving gluons, but only approximately true in the 
four-quark scattering. 

The effect of the Regge excitations can be quantified ^l^IIll through the presence 
of a common form factor in the two-to-two parton scattering amplitudes, which is 
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called the Veneziano form factor and is written in terms of the F-function: 

Vi^,u)^^-^^^lim^^^, (36) 

^ ^ r(i - t/M^) ^ ' 

where s, t, u are the usual Mandelstam variables. The physical content of the 
Veneziano form factor is revealed by an expansion in terms of s-channel poles. 
Each such pole represents a virtual Regge resonance, with mass \/nMs. For the 
purpose of resonances in the dijet spectrum, only the first-level (n = 1) excitation 
is relevant, while the string mass Ms is the only free parameter. The exact values 
of the cross sections depend also on the color factors and spin values of the excited 
states. 



2.2.8. i?6 Diquarks 

In the context of superstring theory in 10 dimensions, anomaly cancellation requires 
that the gauge group is £'§ x £'§. Certain models for the compactification of the 
additional 6 dirnensions, predict that the grand unification symmetry group is iS'g'^- 
The Eq models contain color-triplet scalar diquarks, D and with charges — | 
and I respectively, which couple to the light quarks u, d. 

The interaction Lagrangian between the EQ diquarks and the light quarks is 
given by 

C = \e^,ku^'\{l - lr,)dW^ + ^e,,fei/^i(l + i5)d-^D-^ + h.c, (37) 

where k are color indices, and A, Ac are parameters of the hyper-potential of the 
general Eq mod el. The squared amplitudes for the diquark decays to light quarks 
are given by^^ 

\M{D -^ud)\^ = 24X^ml, {MiD" ud)\^ ^ GXlml^. (38) 
The corresponding widths are: 

Td = axMo, rz5c = ^axJ^D", (39) 
with a\ — A^/47r, ax^ — Ac/47r. 



2.2.9. Color octet technirho 

Technicolor models predict a rich spectrum of technirho vector mesons (p-r) that 
decay predominantly into technipion (ttt) states. In most of the models, pr is either 
heavy, leading to small production cross section, or its decay width to ttt is very 
large, leading to broad peak structures on t he d ijet mass spectrum. However, in the 
context of the "walking technicolor" models ^21!^ the technifermion chiral condensate 
acquires large values, resulting in ttt with large mass, such that the pr — tttttt 
decay is kinematically suppressed. Furthermore, if the technirho is a color mul- 
tiplet, its production cross section in hadron collisions through the qq,gg — ^ pt 
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subprocesses is enhanced. In retrospect, giving the kinematically suppressed decay 
to technipions, px decays almost 100% to partons, appearing as a narrow resonance 
in the dijet mass spectrum. 

It should be noted that the phenomenology of the walking technicolor models 
can be fairly complexly, with one reason being the fact that the px states are mixed 
with ordinary gluons. Because of the matrix form of the corresponding propagator, 
the predicted resonance on the dijet mass spectrum appears higher than the pole 
mass. Also, the resonance cross section cannot be calculated independently of QCD, 
to which it is coupled. For the experimental searches, a simple benchmark model 
is considered ^231^ which predicts a single "bump" on the dijet mass spectrum. The 
corresponding parameters are listed below: 

• standard topcolor-assisted-technicolor (TC2) couplings 

• degenerate technirhos: M(pii) = M(pi2) = M(p2i) = M(p22) = M{p) 
(the pole mass) 

• mixing among the technirhos: A/g — 0, which reduces the mass shift from 
the pole mass, 

• octet technipion mass: M{tt22) = which suppresses the decay px — >■ 

• singlet technipion mass: M(7r22) = Af(7r22)/2, 

• coloron mass: M{V8) — > oo , so that the coloron does not affect the cross 
section, 

• mass parameter Ms — which keeps the branching fraction of the 
process px — > Qt^x small and leads to narrow resonance 



2.2.10. Benchmark Models 

The parton-parton resonance models presented in the previous sections involve lim- 
ited number of free parameters each. The experimental searches traditionally con- 
sider benchmark models, with certain parameter assumptions, which are then used 
to set limits on the masses of the corresponding resonances. Below is a summary of 
the benchmark models: 

• Axigluons: the number of quark flavors to which the axigluon can decay is 
set to = 6, corresponding to the known quarks. 

• Colorons: the number of quark flavors to which the coloron can decay is set 
to TV = 6, and the gauge boson mixing angle is set to cot 9 = 1. 

• Excited quarks: standard model couplings are assumed [fs = f = f = 1) 
and the compositeness scale is set equal to the excited quark mass A = M* . 

• RS graviton: the ratio k/Mpi is set to k/Mpi = 1. 

• W , Z' : standard model couplings are assumed. 

• Eq diquark: electromagnetic coupling constants are assumed aA = ax^ = 

Oie- 
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• Color octet scalars: the gauge coupling is set equal to the QCD coupling 
{k — 1), and the characteristic scale of the interaction is set equal to the 
resonance mass A — M . 

Table[2]sumniarizes the basic properties of the resonances discussed in this review. In 
particular, the decay widths are approximate values for the benchmark models, since 
they also depend on the running of as which should be evaluated at a scale equal 
to the resonance mass. For string resonances the decay width varies significantly, 
depending on the spin and color quantum numbers of the resonances. 



Table 2. Summary of resonances considered in this review. 



Resonance 


Symbol 




Color 
Multiplet 


r/(2»nfl) 


Decay 
Channel 


excited quark 


q* 


1/2+ 


triplet 


0.02 


19 


axigluon 


A 


1+ 


octet 


0.05 




coloron 


C 


1- 


octet 


0.05 


gq 


RS graviton 


G 


2- 


singlet 


0.01 


qq, 99 


-Be diquark 


De 


0+ 


triplet 


0.004 


ud 


color octet scalar 


Ss 


0+ 


octet 


0.04 


99 


color octet technirho 


PT 


1- 


octet 


0.01 


99,99 


heavy W 


w 


1" 


singlet 


0.01 


9iq2 


heavy Z 


Z' 


1- 


singlet 


0.01 


99 


string 


s 


various 


various 


0.003 - 0.037 


99,qg,99 



2.2.11. Model Calculations 

The exact LO calculations of the cross sections and the decay widths of the var- 
ious resonances involve all the relevant Feynman diagrams associated with each 
Lagrangian. In practice, the experimental searches presented in this review are fo- 
cused on narrow resonances, which appear as "bumps" on a steeply falling dijet 
mass spectrum. In all the cases, it is the s — channel decay mode of the resonances 
which produces a " bump" . 

The cross section of a resonance decaying through the s — channel is given by 
the Breit-Wigner expression: 

.(,„) (1 + 2 ^ « ^ 3 + 4) . 16.^ X ^" ; "'^y + ^' , ,40) 

where and F^ are the mass and the total width of the resonance, respectively, 
r(l + 2 ^ R) and T{R 3 + 4) are the partial widths for the creation and the 
decay of the resonance to the specific final state. The spin and color multiplicity 
factor Af is 

(41) 
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where -/Vs^ , Ns-^ 2 the spin multipUcities of the resonance and the initial state 
particles, while Cr and C1.2 are the corresponding color factors. The cross section 
above arises after integrating over cos9* . The angular dependence of the cross sec- 
tion, in the s — channel decay mode, is determined by the spin of the resonance and 
the spin of the final state particles. It should be noted, that in all resonance cases 
decaying to two partons, the angular dependence is expressed as a polynomial of 
cos 6**, as opposed to the dominant QCD background, which exhibits a. t — channel 
pole at cos 9* — 1. More specifically, the angular distributions of the various reso- 
nances, in the s — channel decay mode, are listed below: 

• Eq diquark, color octet scalars: F(cos9*) = const. 

• excited quark: F{cos6*) ^ 1 +cos9*, which becomes F(| cos 6**1) = const, 
(odd in cos 9*). 

• axigluon, coloron, W, Z' : F{cos9*) 1 + cos'^ 9* . 

• RS gravitons: F{gg — > G — > qq) = F{qq — > G — > gg) ^ 1 — cos'*0*, 
F{gg G ^ gg) ^ 1 + Gcos^ei* + cos^ 9* , and F{qq G ^ qq) ^ 
1 - 3cos2 6i* +4cos4 9*. 

where F(cos6'*) = da /d cos 9*. 

In practice, experimental searches impose kinematic constraints on the scattering 
angle 9*, such that the QCD background is suppressed. In this case, the Breit- 
Wigner partonic cross section is written as: 

a{m) = 2 ' (42) 

(m^ — ) -I- mj^T'j^ 

where BR is the branching fraction of the subprocess, and .4cose* is the acceptance 
after the cos9* cut. If the resonance is sufficiently narrow (Tfj << rrifl,), the narrow- 
width approximation holds: 

— — « ^^Sim' - ml). (43) 

Finally, the hadronic cross section in the narrow-width approximation is derived; 
<yhad{mR) = 167r^ X TV X Aose* x BR x 

where the parton luminosity ^ is calculated at t — ml/s, and constrained in the 
kinematic range [yrmn^Vmax]- 

3. Experiment 

In this review paper we consider all searches that used the dijet mass or px spec- 
tra to search for dijet resonances. The searches considered are listed in table |3] in 
chronological order, along with a summary of the energy and luminosity of the 
dataset, and the techniques of the search. The searches by the UAl and UA2 exper- 
iments used data from the proton anti-proton collisions at the CERN SppS collider 



ldL{y 
s dr 



X—, (44) 
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at a center-of-mass energy of 0.63 TeV. The searches by the CDF and DO experi- 
ments used data from proton anti-proton coUisions at the Fermilab Tevatron at a 
center-of-mass energy of 1.8 and 1.96 TeV. The searches by the ATLAS and CMS 
experiments used data from proton proton coUisions at the CERN Large Hadron 
CoUider at a center-of-mass energy of 7 TeV. The search techniques hsted in table |3] 
are discussed in the fohowing sections. 

Table 3. Searches for dijet resonances at hadron colliders. For each search we list the experiment, 
year of data publication, center-of-mass energy, integrated luminosity, techniques used to define the 
resonance shape and the background, dijet mass range of the data, the cut applied on the center of 
mass scattering angle, and the primary reference for the search. 



Expt. 


Yr. 




f L dt 


Resonance 


Background 




cos 6* 


Ref. 






(TeV) 


(p6-i) 


Shape 


Shape 


(TeV) 


Cut 


# 


UAl 


86 


0.63 


0.26 


BW e Gaussian 


LO QCD 


.07 - 


0.3 




m 


UAl 


88 


0.63 


0.49 


BW e Gaussian 


LO QCD 


.11 - 


0.3 


bins 




CDF 


90 


1.8 


0.026 


BW © Gaussian 


LO QCD 


.06 - 


0.5 




M 


UA2 


90 


0.63 


4.7 


Gaussian 


Fit Func. 


.05 - 


0.3 




EH 


CDF 


93 


1.8 


4.2 


BW © Resolution 


LO QCD 


.14- 


1.0 




M 


UA2 


93 


0.63 


11 


Gaussian 


Fit Func. 


.05 - 


0.3 


.60 


M 


CDF 


95 


1.8 


19 


Pythia © 


Sim. 


Fit Func. 


.15 - 


0.9 


.67 




CDF 


97 


1.8 


106 


Pythia © 


Sim. 


Fit Func. 


.18 - 


1.0 


.67 


EH 


DO 


04 


1.8 


109 


Pythia © 


Sim. 


NLO QCD 


.18 - 


1.2 


.67 


M\ 


CDF 


09 


1.96 


1130 


Pythia © 


Sim. 


Fit Func. 


.18- 


1.3 




M 


ATLAS 


10 


7 


0.32 


Pythia © 


Sim. 


Fit Func. 


.20 - 


1.7 


.57 


M 


CMS 


10 


7 


2.9 


Pythia © 


Sim. 


Fit Func. 


.22 - 


2.1 


.57 


M\ 


ATLAS 


llw 


7 


36 


Pythia © 


Sim. 


Fit Func. 


.50- 


2.8 


.57 




CMS 


11 


7 


1000 


Pythia © 


Sim. 


Fit Func. 


.84- 


3.7 


.57 


EH 


ATLAS 


lis 


7 


1000 


Pythia © 


Sim. 


Fit Func. 


.72 - 


4.1 


.54 





3.1. Resonance Shapes 

Narrow resonances are those whose observed width is dominated by experimental 
resolution, where the natural width of the resonance is small in comparison. Searches 
for narrow resonances are prevalent in high energy physics, because they in principle 
only require the use of one shape for the signal, dominated by the experimental 
resolution. This is an approximation, as we will see explicitly in this section. 

Searches for dijet resonances have modeled the resonance line shape using either 
an analytic or a Monte Carlo technique, as shown in Table |3l 

The analytic technique used by the older searches started with a natural shape 
for the resonance which was then smeared with the experimental resolution func- 
tion. That natural shape was either a Breit-Wigner (BW) for wide resonances, where 
the natural width could not be ignored, or a delta-function for narrow resonances. 
The principle of the technique is the idea that the resolution can be understood 
separately from the natural line shape, and the two can then be convolved. In prac- 
tice things were often done this way in older searches because of limitations in 
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the MC modeUng or the detector simulation. The search using 1986 data from 
UAl 1^ considered the entire spectrum, with QCD and resonances determined ana- 
lytically and summed together, and then smeared with the experimental resolution, 
o-(pt)/pt = 0.53/^^+0.05. The UAl search in 19881^ smeared a BW shape in 
a toy MC with the Gaussian experiment al di jet mass resolution, cr/m = 0.11, and 
checked the process with the Isajet MC and a full UAl simulation. The UA2 
search in 1990 ^23, and again in 1993 considered only narrow resonances and 
used the Gaussian experimental resolution for the resonance shapes. The resolution 
varied from 10.7% at 80 GeV to 8.4% at 300 GeV for W qq resonances, and UA2 
noted that the resolution was worse for a gg resonance, 13.2% at 80 GeV 1211. The 
CDF searches in 1990 "^^ and 1993 coherently summed an axigluon resonance 
shape to the QCD background and then smeared with the experimental dijet mass 
resolution. The resolution was cr/m — 0.68/^/m + 0.065 in 1990, and in 1993 it 
was parameterized in a inore complex way and also included the effects of QCD 
radiation from HerwigSO!. The 1990 search was for narrower axigluons, N = 5 and 
10, and the 1993 search was for wider axigluons, A'^ = 10 and 20. The 1993 paper 
also considered a more generic search for three widths, one narrow (F = 0.02M) 
and two wide (F O.IM and 0.2Af). 




Mjj (GeV) (Qev) 

Fig. 3. Resonance shapes from UA2 in 1993. (left) Z' and (right) excited quark simulations 
at a mass of 200 GeV with Gaussian fits overlaid. Reprinted from Ref. 1291 Copyright 1993, with 
permission from Elsevier. 

A ni ce early illustration of narrow resonance shapes was provided by UA2 in 
1993 l^ni and is reproduced in Fig. [31 The Gaussian shapes used to model the sig- 
nal repr esen t reasonably well the central core of the simulated distribution using 
Pythia 1311^ but there is also a long tail at low mass due to final state radiation. 
UA2 accounted for the tail using a signal efficiency, and stated that "this tail is 
re-absorbed in the parameterization of the QCD background". The tail at low mass 
is generally negligible compared to the QCD background, so a Gaussian was an ad- 
equate approximation. Nevertheless, all subsequent searches used the full resonance 
shape including the radiative tail to fit for resonances and set limits. 

The MC technique used by the more recent searches modeled the resonance 
line shape using Pythia with either a fast or full detector simulation. The CDF 
searches in 1995 and 1997 used an excited quark decaying to qg to model 
the resonance shape for all dijet resonances. Uncertainties on the amount of final 
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Fig. 4. Resonance shapes from DO in 2004 and CDF in 2009. (left) DO simulation of three models 
of a 500 GeV resonance from Ref. 32, Copyright 2004, and (right) CDF simulation of four models 
of an 800 GeV resonance from Ref. |33I Copyright 2009 by the American Physical Society. 

state radiation were included by changing the amount of the "radiation tail" at 
low mass by ±50%. The Ga ussian cor e of the dijet mass resolution was reported as 
approximately 10% at CDF |30 | 31 | 33 |^ j^^j. ^.j^g resolution did depend on dijet mass. 
The DO search in 2004^2!^ and the CDF search in 200923], used independent shapes 
for different models of narrow resonances available in Pythia, as shown in Fig. |4l 
The shapes varied primarily due to the different amounts of final state radiation, 
with the low mass tail and the width usually increasing with the number of gluons 
in the final state. 

The CMS search in 2010^^ introduced independent narrow resonance shapes for 
each generic type of final state: quark-quark, quark-gluon, and gluon-gluon shown 
in Fig. [5l They reported in Ref. 35 that for resonance masses between 0.5 and 2.5 
TeV the Gaussian resolution of the core of the distribution "varies from 8% to 5% 
for qq, 10% to 6% for qg, and 16% to 10% for gg ". These were estimated by fitting 
mainly the peak and high- m ass edge to a Gaussian To reduce the radiation 
tail the CMS search in 2011 introduced wide-jets with an effective radius of 
Ai? =1.1. Wide-jets also improved the Gaussian component of the resolution which 
in 2011 was a/M = 1.31/v^ + 0.018 for qq, a/M = 1.56/\/M -h 0.027 for qq, and 
(7 /M — 2.09/v'M -f 0.015 for gg resonances, with resonance mass M given in GeV, 
and between 500 and 3500 GeV The improvement is largest for gluon-gluon 
resonances, as shown in Fig.[5l where the shape using wid e-jets is compared with the 
shape using narrower jets from the anti-fc^ algorithm'^ with a distance parameter 
i? = 0.7. 

The ATLAS search in 2010 considered the explicit dijet resonance shape of 
the excited quark model from Pythia, "including all possible SM final states, which 
were dominantly qg but also qW , qZ , and qj". The reported Gaussian resolution for 
an excited quark ranged from 11% to 7% for q* masses between 0.3 and 1.7 TeV^SS, 
and the reported detector level resolution a„i../mjj was "5% at 1 TeV, 4.5% at 
2 Te V, and asymptotically approached 4% at masses of 5 TeV and above" ESI. In 
2011 ESISSI ATLAS reused the q* shapes to set limits on axigluons, after noting 
that "the axigluon and q* signal templates result in very similar limits ". In summer 

OQ 

2011 ATLAS introduced a shape for gg resonances. The ATLAS searches in 



Searches for Dijet Resonances at Hadron Colliders 19 



20lll '>o i 'J° l also considered generic limits for Gaussian shaped resonances of varying 
width, not just narrow resonances. 

The shapes of narrow resonances presented in the literature are all well described 
by a Gaussian core and a radiation tail, but the different resonance shapes in a 
single experiment also demonstrates that the concept of narrow resonances is often 
approximate. 
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Fig. 5. Resonance shapes from CMS. (left) Simulations of 1.2 TeV resonances shown for the 
different possible final state partons, from Ref. 1351 Copyright 2010 by the American Physical 
Society, (right) Simulation of 2 TeV resonances reconstructed with two different types of jets, from 
Ref. 1371 Copyright 2011, with permission from Elsevier. 



In addition to radiation, another parton-level physical mechanism produces a 
significant tail at low mass, further perturbing the narrow resonance search tech- 
nique: the parton luminosities at low mass are larger than at the resonance pole 
mass, and significantly lift the natural Breit-Wigner tail of the resonances, in the 
cases where the PDFs are falling off rapidly with increasing dijet mass. This fre- 
quently happens at resonance masses approaching the kinematic bound for both sea 
quarks and gluons, and the size of the effect increases with the resonance intrinsic 
width, even for resonances normally considered narrow. The extreme end of this tail 
due to the PDFs is sometimes suppressed in the searches by requiring the partons 
to be have mass close to the pole mass, within a few standard deviations on the dijet 
mass resolution. This is generally a reasonable solution for the shapes, as the QCD 
background overwhelms the signal at low dijet mass. However, the way that this 
tail from PDFs is handled can significantly affect the total resonance cross section 
quoted for specific models, as we discuss in [Appendix A[ 

The narrow resonance search technique, where a single resonance shape domi- 
nated by experimental resolution is used to model the effect of all narrow resonances. 
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is only strictly applicable when the half- width of the resonance, r/2 is significantly 
less than the experimental Gaussian resolution a. The model half-widths are listed 
in tabled and vary from about 0.4% to as large as ag/l for axigluons and colorons. 
For the searches at UAl, UA2, CDF, and DO, as/2 is about a factor of 2 less than 
the dijet mass resolution of roughly 10%, and the half-widths of the other models 
are significantly less than 10%, so the narrow resonance technique should be appli- 
cable. However, the searches at CMS and ATLAS are now looking for very massive 
resonances, producing very energetic jets measured with better resolution. For an 
axigluon or coloron of mass 3 TeV the half- width is as (3 TeV)/2 = 3.9%. For a res- 
onance mass of 3 TeV, the Gaussian dijet mass resolution at CMS in the qq channel 
in 2011 is 4.2%, and at ATLAS at 3 TeV the dijet mass resolution quoted above 
is roughly the same. In addition, both CMS and ATLAS have also included a long 
tail to low mass due to radiation, which is a part of the modification of the natural 
line width into an observed line width, so the experimental resolution is somewhat 
worse than 4.2%. Nevertheless, the natural half-width of 3.9% is comparable to the 
experimental resolution, and so the widest resonances ATLAS and CMS have con- 
sidered, like axigluons and colorons, are beginning to push the boundaries of the 
narrow resonance classification. 



3.2. Dijet Data and QCD Background 

The heart of the search for dijet resonances is the measurement of the dijet mass 
distribution and the estimation of the background. Unlike many other searches in 
high energy physics, the search for dijet resonances is completely dominated by a 
single background. The observed dijet mass distribution comes from the dominant 
process in hadronic collisions: 2 -> 2 scattering of partons predicted by perturbative 
QCD. 



3.2.1. Angular Requirement 

The event selection requirements of each search can only be understood by exam- 
ining the dijet production in QCD. Most experiments chose a fiducial region in the 
experiment to measure dijet production, limiting the pseudorapidity, 77, of each jet 
to a central region. In addition, as shown in table [3l many of the searches place 
a cut on the center of momentum frame scattering angle | cos 9* \ Both of these 
selections, but the | cos 9* \ one in particular, are designed to enhance the dijet res- 
onance signal and suppress the QCD background. This is because QCD production 
of dijets at high mass is dominated by i — channel production, with an angular 
distribution that is approximately Rutherford scattering (Section [2.1.4p : 



dcos9* f2 (l-cos0*)2' ^ ' 

''For some experiments we have translated their cut on dijet Ar; or A.y to the equivalent cut on 
cos 6* tanh{Ar)/2), which follows from Eq. [5]and A»7 Ki Ay. 
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peaking strongly at cos 9* = 1. In practice, only |cos6'*| is measured, because the 
partons in the final state emerge as jets and are experimentally indistinguishable. 
In the s — channel mode, dijet resonances lead to angular distributions that are 
much flatter in cosO* than QCD, with the exact angular distribution depending on 
the spin of the resonance and the final state partons (Section l2.2.1ip . The optimal 
I cos 6** I cut is determined by maximizing the integrated signal, S, over the square 
root of the integrated background, ^/B, often by using a Monte Carlo. However, 
one can estimate the optimal cos e* cuiEl, using equation |45] for the background and 
a signal flat in | cos 6* \ , resulting in the optimal cut | cos 0\jj^-^ I = 1/V3 = 0.58. 
Table [3] shows that all experiments which employed a cut on | cos 9* \ used a similar 
value in order to optimize their searches for dijet resonances. 

3.2.2. Background Models 

Searches for dijet resonances have chosen to model the background with either a 
QCD calculation or with a background parameterization, as shown in table [H No 
search has used a Monte Carlo simulation to model the QCD background, even 
if they used a Monte Carlo to model the resonance signal. This is because Monte 
Carlo may not model the spectrum shape accurately enough. QCD calculations are 
also challenged by high statistics jet data. In the cases where leading order (LO) 
QCD was used to model the spectrum shape, the normalization of LO QCD was al- 
ways adjusted to obtain agreement with the data. Despite the considerable progress 
in next-to-leading order (NLO) QCD calculations over the last 20 years, only one 
experiment has successfully used NLO QCD to model the dijet background 
and remarkably was able to use the absolute normalization of the calculation. Even 
when normalization and shape agree between QCD and the data, there are still con- 
siderable theoretical uncertainties (PDFs, renormalization scale, non-perturbative 
effects, ...) and experimental uncertainties (jet energy scale, resolution smearing, ...) 
which can produce significant differences between data and QCD that can poten- 
tially mimic a resonance signal. For these reasons, the majority of searches simply 
used a fit function to parametrize the QCD background. This is a traditional method 
to search for a large resonance signal. It also takes advantage of the fact that the 
QCD background always produces a smooth and monotonically decreasing spec- 
trum with high statistics that can simply be fit. This methodology makes use of 
Occam's razor noting that if the data can be well fit with the background fit 
function alone, then there is no significant evidence for signal, and the data is com- 
patible with background alone. Hence it is simplest to use the data itself for the 
background. Nevertheless, most of these searches are aware that this argument has 
its limitations, for example on the low statistics tail at high dijet mass, and have also 
compared the data with a QCD Monte Carlo just to confirm that dijet mass data 

•^Integrate both signal and background from — cos d*j^fj^x *° ^MAX ^'^'^ then solve 
d{S /\/B)/d cos eij^x = 0. 
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are at least roughly compatible with QCD alone. Experimenters have also n oted 
other potential biases when using a parameterization to model the background SSI. 



3.2.3. Data from the CERN SppS Collider Experiments 




Fig. 6. Data from UAl in 1986 and 1988. (left) Inclusive jet px spectrum using 260 nb~^, 
originally from Ref. I24l in 1986, compared to predictions of QCD and axigluons decaying to TV = 5 
flavors of quarks from Ref. 1121 Copyright 1988 by the American Physical Society, (right) Dijet 
mass spectra using 490 nb~^ in four bins of | cos 9* | compared to predictions of QCD and a vector 
resonance with width F = O.IM from Ref. |25i Copyright 1988, with permission from Elsevier. 



The earliest search for dijet resonances at hadron colliders was aty pical 

in many respects. In Fig. |6]we reproduce the comparison of UAl 1986 datcP^from 
the inclusive jet pr spectrum for \r]\ < 0.7 with a leading order calculation of QCD 
plus the axigluon model of a dijet resonance. This comparison was published by 
theorists 1^ in the paper which introduced the theoretical cross section for axigluon 
production. The analysis was challenging, as the jet px data had been unsmeared 
for the effects of jet resolution by the UAl collaboration, and had to be resmeared 
by the theorists in order to make a valid comparison with the expected smeared 
shapes of axigluon resonances in the UAl detector. Figure |6] also shows a search 
from UAl in 1988 which used the dijet mass spectrum, and again compared 
with LO QCD normalized to the data. Both searches using UAl data required a 
multiplicative factor of 1.5 to normalize LO QCD to the data. The UAl 1988 data 
extended to a dijet mass of 0.3 TeV and agreed with the QCD background. The 
UAl analysis of the dijet mass spectrum was unique in that it analyzed the data in 
four bins of | cos 9* \ instead of simply applying a single cut on | cos 9* \ . 
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Fig. 7. Data from UA2 in 1990. (left) The dijet mass spectrum using 4.7 pb~^ and (right) the dif- 
ference between the data and the fitted background. The fitted signal from W and Z decays to dijets 
is shown by the solid curve. From Ref. 1271 with kind permission from Springer Science+Business 
Media. 



The UA2 collaboration in 1990 published the only observatior0 in hadron colli- 
sions ^23 of an s — channel dijet resonance: the W and Z boson, whose mass was al- 
ready known from its previous discovery in decays to leptons. In Fig. [7] we reproduce 
the dijet mass spectrum before and after the subtraction of the QCD background. 
The data after the subtraction is well fit by a single peak with shape, accounting 
for both the W resonance and the smaller rate for the Z resonance. The QCD back- 
ground came from a parameterization of the differential rate as a function of dijet 
mass, m, 



dm 



^0 ^-ivim+psr 



with four parameters poi Pii P2 and p^. The fit of this parameterization to the data 
in the full dijet mass region 48 < m < 300 GeV had a probability of only 1%, but 
excluding the dijet mass region 70 < to < 100 that contained the W and Z boson, 
and redoing the fit, gave a probability of 78%. The clear peak after subtracting this 
second background fit is quite significant, as shown in Fig. [71 The search for the 
dijet decays of the W and Z benefited from apriori knowledge of the existence and 



"^Latcr observations of W decaying to dijets always involved more complicated processes, such as 
top quark decay. 
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the mass of the resonance. Otherwise, the 1% probabihty of the first background 
fit would not be sufficient disagreement with the background to claim a discovery. 
Nevertheless, this provided an important confirmation of the W and Z resonances in 
the dijet channel. The background parameterization was also used outside the mass 
region of the W and Z to search for higher mass resonances in 19901^, and again in 
1993!^. These two papers led the field by introducing searches with parametrized 
background shapes. It is unfortunate that no figures exist comparing either UA2 
analysis for m > 0.12 TeV to the background parameterizatior^. Both UA2 analyses 
used dijet mass data in the region 0.05 < m < 0.3 TeV to search for resonances, 
with the 1990 search requiring |ry| < 0.7 and the 1993 search requiring | cos 6**1 < 0.6. 

3.2.4. Data from the Fermilab Tevatron Collider Experiments 




Fig. 8. Data from CDF in 1990 and 1993. (left) Dijet mass spectrum using 26 nb~^ compared 
to predictions of QCD and axigluons decaying to A'^ = 5 (upper) and N = 10 (lower) flavors of 
quarks from Ref. 1261 Copyright 1990, and (right) using 4.2 pb-i for Af = 10 (upper) and Af = 20 
(lower) from Ref. |28| Copyright 1993 by the American Physical Society. 

The two earliest CDF searches in 1990 and 1993 are very similar, and 
model the background roughly like the UAl search!^. Dijet mass data with |7/| < 0.7 
from these CDF searches are shown in Fig. [5] compared to predictions from the QCD 
background and axigluon signals. In 1990 CDF searched for dijet resonances, using 
26 nb^i, with data in the mass interval 60 < m < 500 GeV. In 1993 CDF used 4.2 
pb~^ of data to search in the dijet mass interval 140 < m < 1000 GeV. Both CDF 
analyses modeled the background using LO QCD calculations smeared with the 

•^Also, no figure with all the 1993 UA2 dijet mass data was ever published: UA2 in 1993 reused 
the data of the 1990 search and added a similar amount of new data 
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dijet mass resolution and normalized to the data at low dijet mass. Both datasets 
were compatible with the QCD prediction, however, CDF reported in the 1993 
paper that "a small excess of events is observed in the data between 350 and 400 
GeV" to explain a mass gap in its axigluon exclusions for N ~ 20. 
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Fig. 9. Data from CDF in 1995 and 1997. (left) Dijet mass spectrum using 19 pb~^ compared 
to a background fit, and simulations of QCD and excited quark signals from Ref. 1301 Copyright 
1995, and (right) using 106 pb~^ from Ref. 1311 Copyright 1997 by the American Physical Society. 

The CDF searches in 1995 and 199?]^ used similar techniques, and modeled 
the background like the UA2 searches Dijet mass data with \r]\ < 2.0 and 

I cos 0*1 < 2/3 from these CDF searches are shown in Fig. [9l Despite having 5 — 20 
times more luminosity than the previous search published in 1993, the 1995 and 1997 
data also extended to about m = 1 TeV because the | cos0*| < 2/3 cut suppressed 
the QCD background. The 1995 search fit the background with a functional form 
containing only three parameters: 

while the 1997 search used a similar form with four parameters: 

^ = ^ (1 - m/yi + psm^sr , (48) 
am TO^'i 

where ^/s is the pp collision energy. These parameterizations were motivated by 
leading order QCD. The term mimics the mass dependence of the QCD matrix 
elements (Eq. [T7|) . and was borrowed from the UA2 parameterization in Eq. HHl 
The term (1 — m/y^Y^ was introduced by CDF and mimics the mass dependence 
of the parton distributions at an average fractional momentum x = m/yG. With 
increased statistics in 1997 the additional ps term was needed in order to model 
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the data at high dijet mass. In addition to the background fit used to search for 
new physics, CDF also compared their data to a QCD prediction from Pythia, 
with the prediction normahzed to the data, on average in 1995, and normahzed 
to the low mass data in 1997. In the 1997 paper CDF noted that the data was 
"above the QCD simulation at high mas s", s imilar to the excess they reported in 
the inclusive jet cross section at high pT The source of the difference was not 
understood at the time, but it has since been attributed to the proton PDFs within 
the QCD prediction. By using a parameterization for the background, CDF was able 
to minimize the effects of considerable experimental and theoretical uncertainties on 
the determination of the background shape. In Fig. [10] CDF showed the difference 
between their data and the background fit, and concluded there was no evidence for 
new physics. In 1995 they quoted the local significance of fluctuations interpreted 
as a resonance mass of 250 GeV (2.3(t), 550 GeV (1.3ct) and 850 GeV (l.Scr), and 
in 1997 pointed out that the single bin near dijet mass 550 GeV was 2.6 a above 
the fit. It was not common practice in those years to estimate the "look elsewhere 
effect" used to quote a global significance of observing local fluctuations by also 
taking into account the probability of observing the fluctuation anywhere in the 
mass range considered. In Fig. [10] CDF also compared the data with the shape and 
normalization of excited quark signals, allowing the excluded mass to roughly be 
estimated by eye directly from the data. 
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Fig. 10. Fit and systematics from CDF in 1995 and 1997. (left) Fractional difference between 
data using 19 pb~^ and a background fit, compared to simulations of excited quark signals. Inset 
shows systematic uncertainty for signal. From Ref. l30l Copyright 1995, and (right) using 106 pb~^ 
from Ref. 1311 Copyright 1997 by the American Physical Society. 

Figure [11] shows DO data from a search published in 2004 The analysis 
required \'q\ < 1 and |cos0*| < 0.67 and measured dijet masses in the range 
0.18 < TO < 1.2 TeV. The DO analysis was unique in using NLO QCD to model the 
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Fig. 11. Data from DO in 2004. (left) The dijet mass spectrum using 109 pb~^ and (right) the 
fractional difference between the same data and NLO QCD from Ref. 1321 Copyright 2004 by the 
American Physical Society. 



background. A QCD calculation from JETRAD ^S. with CTEQ6M HSl PDFs and 
renormalization scale ^ = Q.bE^^^-^ was smeared with the measured dijet mass res- 
olution. It agreed remarkably well with the data. The /DF was 25/25, considering 
only statistical uncertainties and without any change to the normalization of the 
QCD prediction. Considering the experimental and theoretical uncertainties, the 
agreement was impressive. This level of agreement between experiment and QCD 
calculation for the dijet mass distribution, both in shape and in normalization, has 
not been seen by any of the other experiments. 

Dijet mass data shown in Fig. [T^] was published by CDF in 2009 1^3]^ the last 
search from the Tcvatron. To use the same dataset for both a QCD measurement 
and a resonance search only I77I < 1 was required, with no cut on | cos9*\, and the 
mass distribution was measured in the range 0.18 < m < 1.3 TeV. The technique 
for modeling the background is the same as the 1995 and 1997 CDF searches, 
however, the following parameterization was developed for an improved fit to the 
high statistics data over a wide mass range: 

da_ ^ po(l-TO/v/i)Pi 
dm (m/ys)(P2+P3in[m/v^) ■ 

This parameterization gave a good fit to the data with jDF = 16/17, and CDF 
used it as the background in the search for resonances. In the same paper CDF 
also compared the data with a full NLO QCD calculation from fastNLO using 
CTEQ6.1mI^, and with renormalization scale /i = p^^'^, and found x^/DF — 
21/21 after taking into account the systematic uncertainties. However, unlike the DO 
search, the shape agreement between data and NLO QCD at CDF without including 
systematic uncertainties was not sufficient to use NLO QCD for the background in 
a search for resonances. 
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Fig. 12. Data from CDF in 2009. (top) Dijet mass spectrum using 1.13 fb~^ compared to a 
background fit, and simulations of excited quark signals, and (bottom) the fractional difference 
between the same data and fit and (inset) a zoom at low mass. From Ref. [33] Copyright 2009 by 
the American Physical Society. 



3.2.5. Data from the CERN Large Hadron Collider Experiments 

The introduction of the Large Hadron Cohider greatly accelerated the pace of 
dijet resonance searches. Data from three ATLAS searches are shown in Fig. [T3l 
all published within roughly a year. The three searches r equ ired \ri\ < 2.5 and 
I cos6'*| < 0.57 for the search in 2010 ^1 and winter 2011 and \r]\ < 2.8 and 
I cos 6** I < 0.54 in summer 2011 I^Sl xhe measured dijet mass range for the three 
searches was 0.2 < to < 1.7 TeV, 0.5 < m < 2.8 TeV, and 0.72 < m < 4.1 
TeV respectively. The lowest dijet mass used increased with each search because 
the jet trigger threshold used for the search increased. All ATLAS searches used 
the last CDF parameterization in Eq. |49] to model the background. As shown in 
Fig. [131 the ATLAS searches all reported the bin-by-bin statistical significance of 
the difference between the data and the background fit, clearly showing that all 
fits were good and there was no evidence for new physics. The ATLAS searches 
introduced the BumpHunter^^ statistical method for finding the most significant 
upward fluctuation, and for quantifying its complete global probability, including 
the "look elsewhere effect" , of coming from the background. The largest upward 
fiuctuations and the corresponding background probabilities were 0.55 TeV and at 
least 51% in 2010, 1.1 TeV and 39% in winter 2011, and 1.25 TeV and 82% in 
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Fig. 13. Data from ATLAS in 2010 and 2011. (top left) Dijet mass spectrum using 315 nh~^ 
compared to a background fit, and simulations of excited quark signals, and the bin-by-bin signifi- 
cance of the difference between the data and the fit from Ref. 1341 Copyright 2010 by the American 
Physical Society, (top right) Same using 36 pb~^ from from Ref. 1361 (bottom) Same using 1.0 
fb~^, and without simulations, from Ref. 1381 



summer 2011. 

The CMS collaboration has published searches for dijet resonances in 2010 
and in 2011 E3. CMS required \r]\ < 2.5 and |cosr| < 0.57, the same as the 
first two ATLAS searches, and measured the dijet mass spectrum in the interval 
0.22 < m < 2.1 TeV in 2010, and in 0.84 < m < 3.7 TeV in 2011. In Fig. [H CMS 
compared data to both a QCD simulation and to a background fit. The uncertainty 
in the jet energy scale decreased from 10% in the 2010 publication to 2.2% in 2011, 
decreasing the uncertainty in the comparison with QCD as shown in Fig. 1141 As 
previously discussed, in 2011 CMS introduced wide jets to improve the resolution for 
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Fig. 14. Data from CMS in 2010 and 2011. (top left) Dijet mass spectrum using 2.9 pb~^ com- 
pared to a background fit, and simulations of QCD and signals from excited quarks and string 
resonances. The Band shows the jet energy scale uncertainty, (bottom left) Ratio of the same data 
to the fit, and compared to the same signals. Left figures from Ref. 1351 Copyright 2010 by the 
American Physical Society, (top right) Same as top left using 1.0 fb~^ and including the bin-by- 
bin significance of the difference between the data and the fit. (bottom right) Ratio of the data 
and the fit compared to excited quark signals. Right figures from Ref. 1371 Copyright 2011, with 
permission from Elsevier. 



resonances decaying to gluons. Following the CDF and ATLAS searches, CMS used 
Eq. |49] to fit the background in the search. The background parameterization gave 
a good fit: ^IDF = 32/31 in 2010 and ^jDF = 27.5/28 in 2011. In 2011 CMS 
adopted the ATLAS style of displaying the bin-by-bin significance of the difference 
between data and fit, just below the data points, from which one can visually see 
that there were no significant fiuctuations. Unlike CDF and ATLAS, CMS did not 
quote estimates of the significance of fiuctuations in either paper. Following the 
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CDF searches, CMS compared their data with the size and shape of various dijet 
resonance signals in a plot of the data divided by the fit, allowing a direct visual 
estimate of the signal mass values excluded. 



3.3. Limits and Models 

Limits on dijet resonances are often performed in three steps. The first step is setting 
limits on the cross section, which can be purely experimental numbers determined 
using a statistical procedure, and can be fairly independent of any model. The 
second step is to evaluate the theoretical cross section for the model as a function 
of resonance mass, which is a completely model dependent process. The third and 
final step is to compare the experimental upper limit on the cross section, with the 
theoretical cross section of the model, and find what resonance masses, if any, are 
excluded. In Table |4] we summarize the mass intervals excluded for the majority of 
models considered in dijet resonance searches. 



Table 4. Excluded mass intervals for resonance models from searches in the dijet mass spec- 
trum. The model parameters used for these exclusions are presented in Section l2.2.10l See the 
text about each exclusion for any exceptions in the model parameters, discussion of the cal- 
culations, and for the additional exclusions of color octet technirhos by CDF and color octet 
scalars by ATLAS. Continuation rows beginning with a " mark arc present when the search 
excluded multiple mass intervals for a given model. 



Expt. 


Year 


Axigluon 
or Coloron 


Excited 
Quark 


W 


Z' 


Eg 
Diquark 


String 






(TeV) 


(TcV) 


(TeV) 


(TeV) 


(TeV) 


(TeV) 


UAl 


1986 


0.13-0.28 












UAl 


1988 


0.15-0.31 












CDF 


1990 


0.12-0.21 












UA2 


1990 






0.10-0.16 








CDF 


1993 


0.22-0.64 












UA2 


1993 




0.14-0.29 


0.13-0.26 


0.13-0.25 






CDF 


1995 


0.20-0.87 


0.20-0.56 










CDF 


1997 


0.20-0.98 


0.20-0.52 
0.58-0.76 


0.30-0.42 




0.29-0.42 




DO 


2004 




0.20-0.78 


0.30-0.80 


0.40-0.64 






CDF 


2009 


0.26-1.25 


0.26-0.87 


0.28-0.84 


0.32-0.74 


0.29-0.63 


0.26-1.4 


ATLAS 


2010 




0.30-1.26 










CMS 


2010 


0.50-1.17 
1.47-1.52 


0.50-1.58 






0.50-0.58 
0.97-1.08 
1.45-1.60 


0.50-2.50 


ATLAS 


2011w 


0.60-2.10 


0.60-2.15 










CMS 


2011 


1.00-2.47 


1.00-2.49 


1.00-1.51 




1.00-3.52 


1.00-4.00 


ATLAS 


2011s 


0.80-3.32 


0.80-2.99 
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3.3.1. Experimental Systematic Uncertainties 

There are comnion systeniatics faced by all the experiments when setting upper 
limits on the cross section of dijet resonances. Here we give an overview of these 
systematic uncertainties and their relative importance for the majority of the ex- 
periments, before we discuss the limits from each experiment below. 

The largest source of experimental systematic uncertainty on the cross-section 
limits is usually the jet energy scale. This is because the QCD background falls 
steeply with increasing dijet mass, and the amount of background underneath a 
potential dijet resonance, at a fixed dijet resonance mass, depends critically on the 
jet energy scale. For a typical spectrum that falls with a large power of the mass, 
dn/dm ^ where A'^ ^ 10 at high masses, a typical 5% uncertainty in the 

jet energy scale leads to an order of 50% uncertainty in the amount of QCD back- 
ground underneath the dijet resonance, which increases the cross section limit by 
roughly 20% in smooth data samples. When there are upward fluctuations in the 
data sample the uncertainty of the limit at a nearby resonance mass can be much 
larger, allowing for the possibility that the fluctuation is attributed to both a res- 
onance and an uncertainty on the jet energy scale. Nearly all experiments correct 
the jet energy to remove non-linearities and non-uniformities in detector response. 
There are potentially more serious uncertainties if these jet energy corrections are 
not continuous, or contain shape deformities, as they create the possibility of man- 
ufacturing a bump, or hiding one. The experiments work hard to insure that this 
does not happen, but they generally do not include any residual uncertainty in a 
quantitative fashion. 

The next largest source of experimental systematic uncertainty is usually the 
estimation of the QCD background. As discussed in Section 13.2.21 assigning sys- 
tematic uncertainties is challenging when a QCD calculation is used to model the 
background. When a simple background parameterization is fit to the data instead, 
two common methods of evaluating the systematic on the background are: trying 
alternate parameterizations, or assigning the statistical uncertainty on the fit pa- 
rameters as a background shape systematic. Either method will frequently give large 
systematics on the parametrized background in the region of high dijet mass, where 
the small amounts of data are not sufficient to constrain the background fit. 

Most experiments report that the uncertainty on the jet energy resolution does 
not have a significant effect on the limits, and in particular is often significantly 
smaller than the effect of the jet energy scale uncertainty. While this may seem 
surprising at first, it can be understood in the same way that the effect of the 
jet energy scale was quantified. On a steeply falling QCD dijet mass spectrum, 
a percentage shift in the mass will have a much larger effect on the amount of 
background underneath the resonance, compared to roughly the same percentage 
widening of the Gaussian resolution. The counter example would be a constant 
background, flat as a function of dijet mass, where a percentage shift in the jet 
energy scale would have no effect on the amount of background underneath the 
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peak, while the same percentage widening of the jet energy resolution would be 
more important. As a conclusion, the resolution uncertainties are negligible because 
the QCD background is falling steeply. 

The astute reader should note that the arguments above are relevant for the 
limits on dijet resonances, which is what most of the searches have reported, and 
the effect of systematic uncertainties would be different if a signal were observed. 

3.3.2. Limits from the CERN SppS Collider Experiments 

The first published exclusion which used UAl data did not use a statistical 
method, or discuss systematic uncertainties. Instead a purely visual exclusion was 
performed. The left plot of Fig. [5] was presented overlaying the observed jet px 
distribution with that from QCD plus the axigluon model. The paper states "We 
see . . . that axigluons in the range Ma = 125 — 275 GeV are ruled out" and "... 
are in contradiction with the data. " 

The search from UAl in 1988 was the first to exclude cross sections for 
resonances at 95% CL, and used the method of maximum likelihood in a fit of the 
QCD background prediction and the resonance shapes to the data. The UAl jet 
energy scale uncertainty was reported as 9%, however, UAl used the agreement 
between LO QCD and the shape of the UAl data at low dijet mass to constrain the 
jet energy scale uncertainty for this search to less than 6% at 95% CL, and then 
included a possible 6% increase in the jet energy scale into the limit. UAl also noted 
that uncertainties in the background estimate from higher order QCD processes were 
not included and could weaken the reported limits. The upper limits shown in Fig. 1151 
were found for narrow resonances (F < O.IM) and wide resonances (F < QAM), 
with the limits corrected for the acceptance of the decay's angular distribution, for 
the two cases of scalar and vector resonances. UAl then compared the wide vector 
resona nce cross-section upper limit to the model cross section from the previous 
search 22] in Fig . [TSl to exclude axigluons in the mass range 150 < via < 310 GeV. 

In 1990 UA2'23j-gported upper limits on the cross section shown in Fig.lTSl Here 
limits at 90% CL were obtained from a fit to the background parameterization and 
a signal, after subtracting the fitted signal for the W and Z resonances. A systematic 
uncertainty on the cross-section upper limit, as much as 21.5%, came from the dijet 
mass resolution uncertainty, which was the dominant systematic. UA2 reduced the 
jet energy scale uncertainty to give only an 11% uncertainty on the cross-section 
upper limit, by calibrating the jets with the observed W and Z peak. UA2 compared 
the upper limits on the cross section to an 0(0^) calculation of the W resonance 
with standard model couplings to exclude the resonance decays in the mass interval 
101 < M < 158 GeV. 

The search from UA2 in 1993^22] dijet resonance signal and the background 
parameterization in Eq. |46] to the data, obtaining the signal cross section and its 
error. From the best-fit signal cross section and its error UA2 calculated the 90% CL 
upper limit including statistical uncertainties only. The systematics were included in 
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Fig. 15. Limits from UAl in 1988 and UA2 in 1990: (left) Upper limits at 95% CL on dijet 
resonance cross section times branching ratio from UAl using 490 nb~^, for both vector and scalar 
particles of two different widths, compared to an axigluon calculation. From Ref. 1251 Copyright 
1988, with permission from Elsevier, (right) Upper limits at 90% CL on dijet resonance cross 
section times branching ratio from UA2 using 4.7 pb~^, compared to a W' calculation. From 
Ref. I27l with kind permission from Springer Science+Business Media. 




Fig. 16. Limits from UA2 in 1993. (left) Upper limits at 90% CL on dijet resonance cross section 
times branching ratio using 11 pb~^, with the nine curves showing systematic variations in the 
limit, compared to a W' cross section calculation, and (right) the excluded mass region for a W' 
in the mass vs. branching ratio plane, from Ref. , 29, Copyright 1993, with permission from Elsevier. 
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a second step and consisted of the uncertainty on the dijet resonances decay width, 
and the uncertainty on the stabihty of the dijet mass measurement between the two 
data samples used from the 1989 and 1990 running periods. For each quantity, x, 
with uncertainty Ax, they performed the cross-section fit with the three different 
values X, X + Ax and x — Ax. Therefore, for the two systematics considered, there 
were nine independent fits and nine different values of the cross-section upper limit. 
For a W signal. Fig. [1^] shows the cross-section upper limit from the nine fits. 
UA2 selected the largest value among the nine as the upper limit on the cross 
section including systematics, and used that to exclude at 90% CL the W in the 
mass region 130 < M < 261 GeV. Also in Fig. \W\ UA2 presented their excluded 
region in the branching ratio vs. mass plane of the W , where a branching ratio of 1 
represents decays to ud and cs quarks, as expected from a right handed W model, 
Wh. This was the only exclusion of a Wn model in dijet resonance searches, as all 
other searches were for a more standard model like left-handed W . Using similar 
techniques, UA2 also excluded at 90% CL Z' bosons decaying to dijets in the mass 
region 130 < M < 252 GeV, and mass degenerate excited quarks (q*) in the mass 
region 140 < M < 288 GeV. Fig. [18] shows the UA2 q* limits in the couphng vs. 
mass plane. UA2 also presented mass limits for a few other variations of the W, 
Z', and q* models ^221 The cross section of the W and Z' models came from an 
0(0^) NNLO calculation, and the cross section for the q* model came from a LO 
Born level calculation. 

3.3.3. Limits from the Fermilab Tevatron Collider Experiments 

The search from CDF in 1990^211 ^Jid not publish upper limits on the cross section. 
A was calculated for a fit of an axigluon plus LO QCD background to the data, 
and a jet energy scale uncertainty ranging from 5% to 9% was included by allowing 
both edges of each dijet mass bin to vary, and CDF used the smallest resulting 
to exclude axigluons with N = 5 decay channels in the mass range 120 < niA < 210 
GeV at 95% CL. In 1993 CDF i^Sj published a table of upper limits on the cross 
section, the first example of publishing generic upper limits, but did not compare 
them with any models in a figure. The 1993 paper directly excluded axigluons with 
iV = 10 decay channels in the mass range 220 < rn^ < 640 GeV at 95% CL, with 
an unspecified statistical test. The paper also set limits on axigluons with A'^ = 20 
decay channels, and presented limits for two choices of PDFs. Both the 1990 and 
1993 papers from CDF used a coherent sum of the axigluon signal and the LO QCD 
background, normalized to the data at low dijet mass. 

In 1995 ESI and 1997211 CDF published cross-section upper limits on dijet reso- 
nances. The statistical technique was described in the prior CDF search for excited 
quarks in the 7 -I- jet and W^-|-jet channel Using a binned likelihood for each 
value of resonance mass, the best fit to the data of the parametrized background 
plus a floating signal were found. For this background the binned likelihood, L, was 
written as a function of the signal normalization a: 
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Fig. 17. Limits from CDF in 1995 and 1997. (left) Upper limits at 95% CL on dijet resonance 
cross section times branching ratio times acceptance using 19 pb— 1, compared to calculations of 
the same observable in six models, from Ref. 1301 Copyright 1995, and (right) the same using 106 
pb~^ , from Ref. l31l Copyright 1997 by the American Physical Society. 



L^TT ^' ^, where ijl, = anAS) + nAB) . (50) 

Here rii is the number of events observed in the data in dijet mass bin i, ni(S) 
is the number of events predicted for the signal using the resonance shape in the 
same bin, and ni{B) is the number of events predicted for the background in the 
same bin. The background parameterization used in 1995 was Eq. |47] and in 1997 
it was Eq. |48l The 95% quantile of the hkehhood distribution in Eq . [50l was found, 
defining an upper hmit on the signal cross section with statistical uncertainties only. 
The dominant sources of systematic uncertainty considered by CDF were the 5% 
jet energy scale uncertainty, the QCD radiation effects on the the mass resonance 
line shape, and the background parameterization, while other sources were also 
included. For each la- shift in the source of systematic uncertainty, the background 
plus signal fit was repeated, giving a new limit. The change in the limit defined a 
la uncertainty in the signal cross section for each source of systematic, and a total 
systematic was derived by adding these in quadrature, as shown in the insets in 
Fig. [ini The likelihood distribution of Eq. [50] was then convolved with the total 
Gaussian systematic uncertainty in the cross section, and the final limit including 
systematics was defined as the 95% quantile of this smeared likelihood. The cross- 
section upper limits on narrow dijet resonances are shown in Fig. 1171 where the 
effect of the small upward fluctuation in the data near a dijet mass of 550 GeV in 
both publications produced a noticeable bump in the upper limit. The upper limits 
were published in a table, and can be compared to the cross section for any model 
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of dijet resonances decaying to two partons with jryj < 2 and |cos0*| < 2/3. CDF 
explicitly did the comparison with a few model cross sections, shown in Fig. 1171 
and the reader can note that the strong processes, like g*, have much larger cross 
sections than the weak processes like W . From this figure CDF excluded in 1995 
and 1997 the mass intervals listed in table |4] for axigluons, excited quarks, W\ Z' 
and Eq diquarks. CDF also excluded color octet technirhos in the mass interval 
320 <M < 480 GeV in 1995 and 260 < M < 470 GeV in 1997. CDF used its own 
LO calculations of the cross section for all models and included 0{as) X-factors 
for the W and Z' models. The calculations for axigluons, W, Z' , and Eq diquarks 
employed the narrow width approximation discussed in section [2.2. Ill In 1997 CDF 
noted that the cross section for the coloron model was always greater than or equal 
to that for the axigluon model, and so all the 95% CL upper limits for axigluons 
also applied to colorons. For the excited quark model, CDF plotted the exclusion in 
the coupling vs. mass plane in Fig. 1181 comparing with the dijet resonance search 
from UA2 in 1993. Figure [18] also shows prior CDF exclusions of q* decays from 
the 7*7 and Wq channel, which were combined with the limit from the dijet channel 
in 1995 to extend the q* mass limit from 560 GeV with dijets alone to 570 GeV 
including 7(7, Wq and gq decays. 
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Fig. 18. Excited quark excluded regions from CDF in 1995 and 1997 including UA2 from 1993. 
(left) Region of the excited quark coupling vs. mass plan excluded by CDF using 19 pb— 1 in the 
qg channel, combined with exclusions in the qW and 57 channel, and compared with the dijet 
exclusions from UA2 in 1993 and also from LEP, from Ref. 1301 Copyright 1995, and (right) same 
using 106 pb— 1, from Ref. 31 Copyright 1997 by the American Physical Society. 

The search from DO in 2004 ^21 published limits on three models of dijet res- 
onances: W, Z' and q* . The handling of statistical uncertainties in the limit was 
similar to the technique used by CDF in 1995 and 1997, which DO noted as a 
Bayesian technique with a flat prior for the signal DO, unlike CDF, applied a 
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truly Bayesian methodology to the treatment of the systematic uncertainties. The 
systematics considered were the jet energy scale (2%), resolution, efficiency, and 
luminosity, with no uncertainty on the background from the NLO QCD calculation. 
All these nuisance parameters had Gaussian prior probability distributions with 
widths given by their uncertainties. Eq. [50] was multiplied by all priors, and then 
integrated over the nuisance parameters to obtain the posterior probability density 
as a function of the amount of signal, from which the limit was found as the 95% 
quantile. This relatively modern procedure likely resulted in a smaller and more 
correct effect of the experimental systematic uncertainties on the limit than the 
conservative procedure used by CDF in 1995 and 1997. In Fig. [19] the upper limits 
at 95% CL on cross section time branching fraction times acceptance are shown 
separately for each of the three models, and compared to the model cross sections. 
DO obtained the q* , W and Z' LO cross sections from Pythia, and also applied 
NLO correction factors of about 1.3 to the W and Z' . The mass intervals excluded 
for the the three models are listed in table [4] Comparing to the CDF limits from 
1997, which used a similar sized dataset from the same running period, we note 
that DO was able to obtain a Z' mass limit while CDF was not, DO W mass limits 
were significantly better than the corresponding CDF mass limits, and the DO q* 
mass limit was a little better than CDF and filled in the CDF gap in the q* mass 
limit around 550 GeV. Comparisons of mass limits are perilous, as the model cross 
sections presented by the two experiments do not agree. For example, at a resonance 
mass of 700 GeV, where the acceptance of the two experiments was the same, the 
DO W model cross section in Fig. [TO] is around 1.3pb, about a factor of 2.6 times 
the CDF W' model cross section, and equal to the DO q* model cross section at the 
same mass. The upper limits on the cross section are generally better to compare 
than mass limits, particularly in regions where the acceptance is the same, and the 
DO cross-section upper limits at 700 GeV in Fig. [TH] are around 0.8 — 1 pb for all 
three models which is better than the CDF limit of 1.3 pb. 

The search from CDF in 2009 published upper limits on the cross section for 
dijet resonances. Similar to the 1997 search, CDF also applied a Bayesian methodol- 
ogy starting from Eq.[50]for the binned likelihood including statistical uncertainties 
only. Unlike the CDF 1997 search but similar to the DO search, the method for 
incorporating systematic uncertainties was truly Bayesian The systematic un- 
certainties considered were the jet energy scale (2-3%), resolution, and luminosity. 
CDF used the parameterization in Eq. [49] to describe the background. To account 
for systematic uncertainties in the background, CDF utilized a profile likelihood 
method, in which for each value of the signal cross section considered the param- 
eters of the background were found again by maximizing the likelihood. This is in 
contrast to the previous CDF analysis in 1997, where the background parameters 
were held fixed as the signal cross section was varied. CDF used the shapes from 
four models of dijet resonances W' , Z' , RS graviton, and g*, available in Pythia, 
to set upper limits at 95% CL on the cross section and noted that the limits get pro- 
gressively worse as more gluons are included in the final state. CDF then compared 
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Fig. 19. Limits from DO in 2004 and CDF in 2009. (left) Upper limits at 95% CL on dijet resonance 
cross section times branching ratio times acceptance from DO using 109 pb~^, compared with three 
model calculations and using resonance shapes from three models: a) excited quarks, b) W , c) Z' , 
from Ref. 1321 Copyright 2004, and (right) same from CDF using 1.13 fb"'^, compared with seven 
model calculations and using resonance shape from four models: a) W , b) Z', c) RS graviton, and 
d) excited quarks, from Ref. |33| Copyright 2009 by the American Physical Society. 

these upper limits to the model cross sections for these four models, from Pythia, 
and to its own LO calculations of the cross section for color octet technirhos, 
axigluons, colorons, and -Eg diquarks in Fig. ll9l The resulting mass limits are shown 
in Table 2] for most models. The CDF cross-section calculations for axigluons, col- 
orons, and _Bg diquarks used the same technique as in its 1997 search, and CDF 
included the customary NLO k- factor of 1.3 for the W and Z' models. The color 
octet technirho cross section came from Pythia, where the model included the 
QCD background, and which had to be subtracted off in order to obtain the net 
resonance cross section CDF excluded color octet technirhos in the mass inter- 
val 260 < M < 1100 GeV in 2009. No mass limits were set on Randall-Sundrum 
gravitons by this search, or by any other dijet resonance search up to this date. 
CDF presented the upper limits in a table for future use, which allowed the CMS 
experiment to compare the CDF upper limits on the cross section to that expected 
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for a string resonance and determine a CDF mass limit of 1.4 TeV on string reso- 
nances jjj^ 2009 CDF published the most stringent mass limits to date on color 
octet technirhos and the dijet decays of Z' . 

3.3.4. Limits from the CERN Large Hadron Collider Experiments 
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Fig. 20. Limits from ATLAS in 2010, winter 2011, and summer 2011. (top left) Observed and 
expected upper limits at 95% CL on excited quark cross section times acceptance using 315 nb— 1, 
compared to predictions of excited quarks with various tunes and PDFs. Shaded bands are Icr 
and 2(T variations in the expected limit. From Ref. 1341 Copyright 2010 by the American Physical 
Society, (top right) Same using 36 pb~^ and in addition comparing to an axigluon prediction. From 
Ref.lMl (bottom left) Same using 1.0 fb ^ and (bottom right) same, except limits and predictions 
are specifically for a color octet scalar resonance decaying to gg, from Ref. 1381 



ATLAS pubhshed limits from three datasets: in 20101^, winter 2011 1^ and 
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summer 2011 with 0.3, 36, and 1000 pb~^ of integrated luminosity. Limits 
were set using a Bayesian statistical technique with uniform prior for the signal 
cross section, and systematic uncertainties were incorporated in a fully Bayesian 
treatment using Gaussian priors. The dominant source of systematic uncertainty 
was the jet energy scale in all three publications: 6%-9% in 2010, 3.2%-5.7% in 
winter 2011, and less than 4% in summer 2011. The background was determined 
using the last CDF parameterization in Eq. 1491 The systematic uncertainty on the 
background was taken from the statistical uncertainty in the fitted parameters, and 
for the winter 2011 search the uncertainty increased with resonance mass from 3%- 
40%. The luminosity uncertainty was 11%, 11%, and 3.7% for the three searches 
respectively. In all three searches the uncertainty in the jet energy resolution had a 
negligible effect on the limits. 

ATLAS published the first dijet resonance searches that included expected limits 
and their variation at the la and 2a levels. These are limits determined from pseudo- 
experiments generated from the smooth background prediction. Expected limits 
vary smoothly as a function of resonance mass, while observed limits have wiggles 
that reflect statistical fluctuations in the data. In Fig. [201 the ATLAS observed and 
expected limits at 95% CL are compared with the calculations of the model cross 
sections. In 2010 the fluctuation at 0.55 TeV mentioned previously produced a 2a 
fluctuation in the observed limit, their largest upward fluctuation. 

In all three publications ATLAS compared their cross-section upper limits from 
a q* shape to the cross section times acceptance for the q* model from Pythia 
including all decays: dominantly qg, but also including roughly an additional 20% 
"dijet" cross section resulting from qW, qZ, and qj. ATLAS explored variations 
in the q* cross section with different Monte Carlo tunes and associated PDFs, and 
chose to quote as the main resuh the MC09 tune ^ with MRST2007 PDF 
In 2011 ATLAS also cornpared its cross-section upper limits to an axigluon calcu- 
lation using CalcHEP again with MRST2007 PDF. We discuss the ATLAS 
calculations more in [Appendix A[ The excluded mass intervals for q* and axigluons 
are listed in table HI and the expected and observed mass limits are compared in 
table [5l In the summer of 2011 ATLAS set cross-section upper limits on color octet 
scalars, which decay to gg and have a wider line shape than q* or axigluons. ATLAS 
compared with the model cross section from a MADGRAPH^^plus Pythia calcu- 
lation to exclude color octet scalars in the mass interval 1.0 < M < 1.92 TeV when 
the mass limit 1.77 TeV was expected. All observed mass limits from ATLAS are 
larger than their expected mass limits, due to downward fluctuations in the data. 
ATLAS in the summer of 2011 has published the most stringent limits to date on 
axigluons, excited quarks and color octet scalars. 

In 2011 ATLAS introduced generic upper limits for wid e dij et resonances based 
on a Gaussian line shape. In winter 2011 ATLAS published a table of their 95% 
CL upper limits on the cross section for dijet resonances with measured Gaussian 
widths of 3%, 5%, 7%, 10% and 15% of the peak mass. In summer 20111211 ATLAS 



42 Robert M. Harris and Konstantinos Kousouris 



J3 

X 

e 

c 
o 

i 

-I 10"' 

_] 

O 

Ln 

10^ 



1000 2000 3000 4000 

Mass, rrig [GeV] 

Fig. 21. Model independent limits from ATLAS in summer 2011. Observed upper limits at 95% 
CL for different Gaussian widths, from Ref. 1381 

updated all the limits, except they dropped the 3% Gaussian width, and published 
them in a table and also in Fig. [5T] Notice that the resonances in Fig. [5T] are not 
narrow, and the upper limits on the cross section increases as the Gaussian width of 
these resonance increases. In summer 2011 ATLAS published detailed instructions 
for how to obtain mass limits on new particles from these generic upper limits. In 
short, they recommended the production of a Monte Carlo sample for the calculation 
of the acceptance at parton level, and for the smearing of a dijet resonance signal 
with the ATLAS dijet resolution, which should then be truncated within ± 20% of 
the resonance peak before summing to find the total cross section. This is the first 
publi cati on of generic limits on wide resonances since the 1993 publication from 

cdfISH. 

While this review is only covering searches in the dijet mass spectrum, we note 
that in 2011 ATLAS published the only search for dijet resonances at hadron col- 
liders using the dijet angular distribution and used it to set limits on excited 
quarks. The excited quark limits that we quote for ATLAS in table |4] and table [5] 
are from the search in the dijet mass spectrum only, not the search in the dijet 
angular distribution. 

CMS published limits in 201oE5] using 2.9 ph'^ and in 201l'^ using 1 fb-^ As 
for CDF, DO, and ATLAS, the handling of statistical uncertainties in the limits used 
a Bayesian procedure with uniform prior on the cross section. Systematic uncertain- 
ties considered were on the jet energy scale (10% in 2010 and 2.2% in 2011), jet 
energy resolution (10% producing negligile effect), luminosity (11% in 2010, 6% in 
2011), and the background. The background parameterization used was Eq. |49l the 
same as used by CDF in 2009 and ATLAS. The uncertainty on the background was 
determined in 2010 by varying the choice of the parameterization similar to CDF 
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Fig. 22. Limits from CMS in 2010 and 2011. (top left) ) Observed upper limits at 95% CL on dijet 
resonance cross section times branching ratio times acceptance using 2.9 nb~^, shown separately 
for gluon-gluon, quark-gluon, and quark-quark resonances, and compared with predictions from 
various models. From Ref. 1351 Copyright 2010 by the American Physical Society, (top right) Same, 
usmg 1.0 fb-i and (bottom) same for only quark-gluon resonances and including expected limits 
and the Icr and 2a variations in the expected limits, from Ref. 1371 Copyright 2011, with permission 
from Elsevier. 



in 1997, and in 2011 by varying the fit parameters within their statistical uncer- 
tainty, similar to ATLAS. The procedure for incorporating systematic uncertainties 
used by CMS in 2010 was the same as the method employed by CDF in 1997, and 
was described as "an approximate technique, which in our application is generally 
more conservative than a fully Bayesian treatment". In 2010 the systematic uncer- 
tainties increased the cross-section upper limits by 17%-49%. In 2011 CMS used a 
fully Bayesian procedure with Gaussian priors for incorporating systematics. The 
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updated procedure, along with the reduced systematics on the jet energy scale, had 
a much smaller effect on the observed cross-section upper limits than in 2010. 

CMS used the shapes of narrow qq, qg and gg resonances to set generic upper 
limits at 95% CL on the cross section times branching ratio into dijets times accep- 
tance, and published the numbers in tables in both 2010 and 2011. In Fig. [22] the 
generic upper limits are compared to the dijet decays of string resonances, excited 
quarks, axigluons, colorons, Eq diquarks, W , Z' , and Randall-Sundrum gravitons. 
The model cross sections were LO calculations using CTEQ6L parton distributions, 
the same calculation as at CDF in 1997, and included the same "K-factors" for W 
and Z'^^. The calculations for axigluons, W, Z', and Eq diquarks employed the 
narrow width approximation discussed in section 12.2.111 The CMS axigluon calcu- 
lation was very different from the calculation done by ATLAS, and the CMS q* 
calculation was slightly different than ATLAS, as discussed in [Appendix A[ The 
excluded regions for all models are shown in table [4j and the expected and observed 
limits are compared with ATLAS in table [5] The differences for a data sample of 1 
fb^^ are discussed in the next paragraph. In 2011, CMS also measured the la and 
2a variations in the expected limit, and showed these for qg resonances in Fig. 1221 
The largest upward fluctuation in the 2011 data was near 2.5 TeV which caused the 
observed limit to be less than the expected limit for both excited quarks and axiglu- 
ons. In 2011 CMS published the most stringent limits to date on string resonances, 
Eq diquarks, and the dijet decays of W . 



Table 5. Expected and observed mass limits at ATLAS and CMS from searches 
in the dijet mass spectrum. 
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The CMS publication in 2011 ^'^ and the ATLAS publication in summer 201lE£l, 
both using an integrated luminosity of 1 fb~^, reported different upper limits for 
the same model as shown in table [5| The greatest difference lies in the observed 
limits for axigluons (CMS 2.47 TeV, ATLAS 3.32 TeV) but there is also a signif- 
icant difference for excited quarks (CMS 2.49 TeV, ATLAS 2.99 TeV). As noted 
above, the observed limits by ATLAS were greater than their expected limits due 
to a downward fluctuation in the data, and the observed limits by CMS were less 
than the expected limits due to an upward fluctuation in the data, so these fluctu- 
ations clearly contributed to the difference in the observed limits between the two 
experiments. The fairest comparison should be in the expected limits. However, 
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there is still a significant difference in the expected mass limits for axigluons (CMS 
2.66 TeV, ATLAS 3.07 TeV) and a smaller but non-negligible difference in the ex- 
pected mass limits for excited quarks (CMS 2.68 TeV, ATLAS 2.81 TeV). For both 
models, the vast majority of the difference in the mass limits appears to result in 
differences in the calculated cross section of the model, which we discuss in detail 
in [Appendix A[ 

The most direct and natural comparison between the performance of ATLAS and 
CMS is in their expected limits on the cross section, as opposed to their expected 
mass limits on models. At a resonance mass of 3 TeV, the ATLAS expected upper 
limit on cross section time acceptance for excited quarks is 0.023 pb, while the CMS 
expected upper limit on cross section times acceptance for qg resonances is 0.028 
pb. For excited quark signals at high masses the CMS acceptance is 14% greater 
than the ATLAS acceptance. After correcting for acceptance the ATLAS expected 
limit on the cross section is about 10% more stringent than the CMS limit, which 
would give an expected mass limit about 1% more stringent. Compared this way 
the ATLAS and CMS performance is similar. This is expected given the comparable 
capabilities of the two experiments with the same integrated luminosity at the same 
center-of-mass energy. 

4. Summary 

Dijet resonance searches at hadron colliders have constrained a rich variety of models 
of new physics. Limits have been set on models motivated by grand unification, 
string theory, technicolor, compositeness, and ideas for new color interactions. Mass 
limits on the majority of models constrained are summarized in table ID Limits 
on many colored resonances, like axigluons, colorons, and excited quarks, are now 
around 3 TeV and limits as large as 4 TeV have been set. In addition to model 
specific limits, the searches have provided model independent limits on the cross 
section for dijet resonances that can be used to constrain future models of new 
particle production. 

The most important factors influencing the sensitivity of dijet resonance searches 
are the center-of-mass energy of the collider and the integrated luminosity of the 
dataset. The largest experimental uncertainty, the jet energy scale, now contributes 
little to upper limits on the cross section. Nevertheless, the experiments must re- 
main vigilant that uncertainties in the jet energy scale do not manufacture a dijet 
resonance signal or hide one. Uncertainties in jet resolution at the experiments have 
usually had negligible effect. The searches are always dominated by statistical un- 
certainties, and the experiments have therefore for the most part done the most 
important thing, which is to rapidly search and publish when the energy or inte- 
grated luminosity have increased significantly. Keeping the search simple and model 
independent has made this easier. 

In the course of time, the techniques and ideas employed by the experiments 
have varied significantly. However, an evolution can be observed towards standard 
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practices, and it is worth noting the most important ones which may be useful for 
future searches as weh. 

The existence of a permanent record of fluctuations in the data for future refer- 
ence, is helpful when experiments publish detailed comparisons of their dijet mass 
data to the background prediction. This was attempted by CDF and DO with a ratio 
plot, but the plot of the bin- by-bin significance of the difference between data and 
background, introduced by the ATLAS collaboration and adopted by CMS, is ideal 
for this purpose. Estimations of the global significance of fiuctuations introduced by 
ATLAS are instructive, and it would be even more helpful if they were accompanied 
by the undiluted estimates of the local significance of the fluctuation. 

Both the ATLAS and CMS experiments have adopted a parametrization from 
CDF to model their QCD background. This is because a parameterization fit to the 
data generally gives a better model of the background than QCD calculations. 

The experimental practice of reporting cross-section upper limits on dijet reso- 
nances, in addition to mass limits on specific models, has been important to allow 
continued use of the data to constrain models (see e.g. Ref 13). We have also shown 
that the cross section limits are important to understand and compare the results 
of the experiments, because differences in the calculation of the model cross sections 
often affect the mass limits significantly, as we saw in the comparison of the CDF 
and DO limits with 100 pb"^ and the ATLAS and CMS limits with 1 Ad'^ Clearly 
it would also be helpful if the experiments used a common method of calculating 
the model cross sections. 

Finally, to understand the results of the experiments and to compare them, it 
has been important to have expected upper limits accompanied by their statistical 
variations, as well as observed upper limits. This practice was introduced by ATLAS 
and adopted by CMS, and is a visible benefit of the modern statistical techniques 
that are now commonly employed in the searches. 

We look forward to future searches for dijet resonances at hadron colliders and 
their discovery of new physics beyond the standard model. 
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Appendix A. Axigluon and Excited Quark Calculations from 
ATLAS and CMS 



Axigluon Calculations from ATLAS and CMS 
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Fig. 23. Ratio of the ATLAS and CMS calculations of the axigluon cross section times branching 
ratio times acceptance. 



Fig. [23] shows the ratio ATLAS / CMS for the calculated cross section times 
acceptance for axigluons decaying to dijets in each experiment. The solid curve in 
Fig.[23lshows the ratio of the cross section times acceptance curves published by the 
two experiments in Fig.[20land Fig.[22]for 1 fb^^. At an axigluon mass of 1 TeV the 
calculations agree, and the ratio is the relative acceptance of the two experiments, 
but the ratio increases rapidly with mass. We will use a 3 TeV axigluon resonance 
as an example, where the ratio of the cross section times acceptance published by 
ATLAS and CMS is a factor of 2.8. The Axigluon models employed by the two 
experiments are the same but there are differences in the way the two experiments 
perform the cross section calculations. The major sources of difference are 

(i) Resonance tail. CMS uses the narrow width approximation which accounts 
for only the cross section at the resonant pole, while ATLAS uses CalcHEP 
and integrates the cross section over a dijet mass window within ±30% of the 
resonance pole (0.7M < m < 1.3M) picking up a large tail at masses beneath 
the pole. The affect of ATLAS integration over the mass window is shown by 
both a CalcHEP calculation (dashed curve in Fig. [23]) and by our own 
calculation of the affect of integrating a Breit-Wigner on the CMS calculation 
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(points in Fig. 1^51) . After taking out differences due to acceptance, this increases 
the ATLAS cross section by a factor of 1.92 for a 3 TeV axigluon according to 
the CalcHEP calculation. 

(ii) Parton distributions. CMS uses CTEQ6L and ATLAS uses MRST07. Compar- 
ing the dashed and dotted curves in Fig. [23l which come from a CalcHEP 
calculation we find this increases the ATLAS cross section by a factor of 
1.37 for a 3 TeV axigluon. 

(iii) Processes included. CMS includes only the process where axigluons are pro- 
duced resonantly in the s-channel, while ATLAS also includes processes where 
axigluons are exchanged in the t-channel which increases the ATLAS cross 
section by a factor of 1.18 for a 3 TeV axigluon. 

(iv) Acceptance. The acceptance of the two experiments for Axigluons is dominated 
by the | cos6'*| cut which is 0.57 for CMS and 0.54 for ATLAS. The combined 
affect of this cut and the ATLAS removal of a narrow region of the detector 
decreases the ATLAS acceptance by a factor of 0.89. 

The differences in the axigluon cross section calculations at ATLAS and CMS 
are the largest source of the difference in the expected limits and requires further 
discussion. The resonance tail at low mass is the largest single effect. How to handle 
its contribution to the new particle cross section, is a common problem faced by 
every resonance search in hadron colliders. It is caused by the decrease of the parton 
distributions of the proton with increasing fractional momentum x ~ m/-y/s, and is 
a particularly severe problem for qq resonances at pp colliders. The narrow width 
approximation used by CMS, and discussed in section [2.2.11[ is commonly employed 
by theorists to calculate the cross section for new particles. It replaces the true shape 
of the resonance with a delta function at the resonance pole, and therefore matches 
well the experimental approximation of both the ATLAS and CMS experiments 
when employing narrow resonance shapes to search for new physics and set limits on 
the cross section. However, the narrow width approximation clearly underestimates 
the true axigluon cross section, a significant fraction of which is off the pole at lower 
masses. The ATLAS choice to integrate over the actual line shape gives a total 
axigluon cross section that may be more realistic. However, this may not match 
as well the assumption of a narrow resonance shape peaked at the pole which was 
used to set cross section limits. The ATLAS choice to include t-channel processes 
in the total cross section again gives a result that may be more realistic than the 
CMS choice to only include s-channel processes. However, it is unclear whether 
the ATLAS and CMS methodology of fitting the data for the background shape 
and normalization effectively absorbs some of the cross section of these additional 
processes into the background, setting limits on only the bump like component of the 
resonances. The same can be said for the remainder of the axigluon tail at masses 
outside the search window which is ignored by both experiments. Further, both 
experiments ignore coherence between QCD and axigluons, which may significantly 
affect the tails of the distribution, and is likely an issue when the tail has a larger 
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cross section than the peak. Finally we note that the different choice of PDFs by 
the two experiment is a source of variation in the mass limits that is easy to remove 
in the future. 

A similar comparison between the model cross sections for excited quarks has 
been done at CMS and ATLAS. Here the calculated values of {a x A) at ATLAS 
at a mass of 3 TeV are only 16% greater than the reported values of {a x B x A) 
at CMS. There are three factors making up this difference: 

(i) Parton distributions. CMS uses CTEQ6L and ATLAS uses MRST07, making 
the ATLAS cross section greater by about 10% at M=3 TeV. 

(ii) Branching Fraction. CMS includes only q* — )• qg decays in the dijet cross 
section, which make up 84% of the branching fraction. ATLAS also includes 
the qW, qZ and 97 decays. This choice increases the ATLAS cross section by 
a factor of roughly 1/0.84 = 20% 

(iii) Acceptance. The acceptance of the two experiments for excited quarks, which 
at resonance mass of 3 TeV is 0.57 for CMS and 0.50 for ATLAS. The ATLAS 
acceptance is smaller by 12%. 
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